Could light harvesting complexes exhibit non-classical effects at room temperature? 



Mark M. Wilde and James M. McCracken 
Electronic Systems Division, Science Applications International Corporation, 
4001 North Fairfax Drive, Arlington, Virginia, USA 22203 



o 
o 

> 

o 
in 



a,: 
cr 



> 
o 



ON 

o 



Ari Mizel 

Laboratory for Physical Sciences, 8050 Greenmead Drive, College Park, Maryland, USA 20740 

(Dated: November 5, 2009) 

Mounting experimental and theoretical evidence suggests that coherent quantum effects play a 
role in the efficient transfer of an excitation from a chlorosome antenna to a reaction center in 
the Fenna-Matthews-Olson protein complex. However, it is conceivable that a satisfying alternate 
interpretation of the results is possible in terms of a classical theory. To address this possibility, we 
consider a class of classical theories satisfying the minimal postulates of macrorealism and frame 
Leggett-Garg-type tests that could rule them out. Our numerical simulations indicate that even 
in the presence of decoherence, several tests could exhibit the required violations of the Leggett- 
Garg inequality. Remarkably, some violations persist even at room temperature for our decoherence 
model. 
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I. INTRODUCTION 



Does quantum mechanics play a nontrivial role in 
evolutionary or molecular processes in biological sys- 
tems? How could this result be true, when biological 
systems interact in environments that are "hot" and 
"wet?" Furthermore, how can we frame the first ques- 
tion rigorously? These questions have been the sub- 
ject of intense debate p], Q and a surge of recent in- 
terest from both th e q uantum information community 
0, 1, i, i, 0, 1, i, El EI El 0, 0, El Eland the exper- 
imental quantum chemistry community [17l EH, El, H3] • 
This active area is sometimes described as "quantum 
biology" [1 EH- Some have suggested that, in addition 
to the possibility of remarkable new insights into biol- 
ogy, there is the potential to harness naturally occurring 
coherent quantum effects in biology to engineer new de- 
vices. 

Two biological processes are playing a particularly im- 
portant role in fueling interest in this subject: energetic 
excitation transfer in the Fenna-Matthews-Olson (FMO) 
li g ht-harvestin g p rotein complex of g reen bacteria ujl \4 , 

El El, El El El El, Hi Eil, HI Ha SI HI M Mm\ 

and the radical ion-pair mechanism in the so-called avian 
compass [1 H H El E3]. For each system, there are 
claims of quantum mechanical behavior in the scientific 
literature. On the theory side, several studies in the 
literature have proposed an open quantum systems ap- 
proach to model the dynamics of these biological pro- 
cesses d |3, El El] and have lent credence to the propo- 
sition that coherent quantum effects play a role in their 
functionality. 

Agreement between quantum theoretical models and 
experiment, however, does not irrevocably demonstrate 
the presence of quantum effects. This point is subtle, 



but it is a logical fallacy, called the "affirmation of the 
consequent," to conclude once and for all that the quan- 
tum biological models are the correct models simply be- 
cause they coincide with the observations of some experi- 
ments (Leggett stresses this point in several of his papers 
[H[3(|). It would be useful to irrevocably exclude cer- 
tain classes of classical models that might apply to these 
biological systems by considering fundamental tests of 
non-classicality. 

One might wonder if it is well-motivated to conduct 
a test for non-classicality on a system such as the FMO 
complex that is clearly microscopic. In fact, is it not 
obvious that the FMO complex should behave accord- 
ing to quantum-mechanical laws, given that the length 
scales and time scales are those to which quantum the- 
ory applies? We argue that applying a test for non- 
classicality is valuable here because a significant amount 
of environmental noise acts on the chromophores of the 
FMO complex. These decohering effects may "wash out" 
quantum-mechanical behavior and make it appear as if 
the FMO complex behaves according to a class of macro- 
realistic models (after all, some argue that classical be- 
havior in general arises due to such decoherence [3l|). 
In some cases (in a high temperature limit), researchers 
actually have modeled the dynamics of the FMO com- 
plex with a classical incoherent hopping (Forster) model 
[H [H H3, Recent results B S, El El suggest 

that a classical incoherent hopping model is insufficient 
to explain the ultra-efficient transfer of energy in photo- 
synthesis, but a test for non-classicality could irrevocably 
exclude the whole class of models to which the classical 
incoherent hopping model belongs. 

In this paper, we frame a test for non-classicality in 
the FMO protein complex, the molecular complex re- 
sponsible for the transfer of energetic excitations in a 
photosynthetic reaction. Quantum chemists have deter- 
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mined a tight-binding Hamiltonian for this simple sys- 
tem [24j. Phenomenological modifications to the stan- 
dard Schrodinger equation have allowed an open quan- 
tum systems model of its dynamics [3, 0, 

mm. M ° ti_ 

vated by experimental results, several theoretical studies 
have computed the efficiency of energy transfer from the 
chlorosome antenna to the reaction center in green bacte- 
ria and have asserted that coherent quantum effects play 
a role. To assess this claim, we frame a test for viola- 
tion of macrorealism, as quantified by the Leggett-Garg 
inequality 34] . In this sense, we are following the pro- 
gram of Leggett outlined in Ref. [29[ and the suggestion 
of Zeilinger in Ref. [l[ that it would be useful to sub- 
ject biomolecules to tests of non-classicality. The Bell 
}35L l36l ] and Leggett-Garg [f33| tests of non-classicality 
are benchmarks that determine whether a given dynam- 
ical system has stronger-than-classical spatial or tempo- 
ral correlations, respectively. Each of these tests provides 
an inequality that bounds the spatial or temporal corre- 
lations of a given system — a violation of the inequality 
implies that the system in question does not behave in 
a classical manner. The advantage of the Leggett-Garg 
test over a Bell test is that a Leggett-Garg test applies to 
a single quantum system, easing experimental difficulty. 
The Leggett-Garg inequality was originally applied to su- 
perconducting quantum systems. It was later considered 
for other systems such as quantum dots [13, [H, [39[ and 
photons [401 ] but, to our knowledge, this is the first appli- 
cation of the Leggett-Garg theory to a biological system. 

We organize the rest of this paper as follows. We intro- 
duce the Leggett-Garg inequality in the next section and 
review a simple example of a two-level quantum system 
that violates it. In Section IlIIl we review the quantum 
dynamical model of the FMO complex (we specifically 
employ the nine- level model of Refs. [!, 0], which is a 
modification of the seven-level model of Rebentrost et al. 
m Refs. [IlGilll). We then discuss several examples 
of dichotomic observables for the FMO complex to which 
we can apply the Leggett-Garg inequality. Section [TV] 
begins our study of the Leggett-Garg inequality and the 
FMO complex. We first study the Leggett-Garg inequal- 
ity with purely coherent dynamics and are able to derive 
analytical results. These analytical formulas allow us to 
determine exactly when coherent dynamics give a viola- 
tion of the inequality. We then exploit these analytical 
results in our numerical simulations of the Leggett-Garg 
inequality and the FMO complex, where we show that it 
is still possible to violate the inequality even in the pres- 
ence of noise and potentially even at room temperature. 
We conclude with some open questions for further study. 



II. THE LEGGETT-GARG INEQUALITY 

The Leggett-Garg inequality applies to any system 
that obeys the postulates of macrorealism. The postu- 
lates of macrorealism for a two-level system are as follows 
00: 



1. Macrorealism per se: A macroscopic object is in 
one of two definite states at any given time. 

2. Noninvasive measurement: It is possible in princi- 
ple to determine the state of the system without 
affecting it or any subsequent dynamics. 

3. Induction: The properties of ensembles are deter- 
mined exclusively by initial conditions (and, in par- 
ticular, not by final conditions). 

It is reasonable to assume that a classical system, in 
principle, should obey the postulates of a macrorealistic 
theory. 

The Leggett-Garg inequality bounds the two-time cor- 
relation functions of three dichotomic observables Q (ti), 
Q{t2), and Q (£3) measured at respective times ti, t 2 , 
and £3. The observables Q (ti), Q (t 2 ), and Q (£3) could 
be the spin of a particle or the location of the trapped 
magnetic flux in a double-well potential as in the original 
application of Leggett and Garg [34| . Let Cij denote the 
following two-time correlation function: 

C id = (Q (<0Q (<,-)>• 

The Leggett-Garg quantity K is the following combina- 
tion of three two-time correlations and a constant: 



K = Ci 2 + C 2 , 3 + C13 + 1. 



(1) 



Note that we can obtain the alternate Leggett-Garg 
quantities 



— Cl,2 + C2,3 — Ci,3 + 1, 

— Cl,2 — C2,3 + Ci,3 + 1, 
C12 — C2,3 — Cl,3 + 1; 



(2) 
(3) 
(4) 



merely by flipping the sign of the respective observables 
Q (*i)i Q (h), and Q (£3). The following Leggett-Garg in- 
equality bounds the Leggett-Garg quantity K when the 
system in question is a macrorealistic system being mea- 
sured noninvasively: 



K > 0. 



(5) 



The last correlation function is to be obtained exper- 
imentally by measuring at times t\ and £3 but refraining 
from measuring at time £2. By comparing Cj.,3 to corre- 
lation functions C\ 2 an d C 2 ,3 obtained in the presence 
of a measurement at t 2 , the Leggett-Garg inequality is 
sensitive to invasiveness in the £ 2 measurement. 

An example of a system that violates the Leggett-Garg 
inequality is a spin-1/2 particle with system Hamiltonian 
H = ijjux /2 and with the observables Q{t\) — Q (£2) = 
Q (£3) = &z- This choice leads to the following value of 
the Leggett-Garg quantity A" in |T]): 

K = cos (2a; At) + 2 cos (wAt) + 1, 

where the parameter At is the uniform time interval 
between the successive measurements of the observable 
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Q. Observe that choosing the interval At = 3it/4uj sets 
K = — \/2+l and leads to a violation of the Leggett-Garg 
inequality in {5). Thus, this quantum system does not 
obey the postulates of a macrorealistic theory when we 
choose the measurement time intervals as given above. 
This violation is perhaps not surprising because a spin- 
1/2 system is a genuine quantum system and "cannot 
have the objective properties tentatively attributed to 
macroscopic objects prior to and independent of mea- 
surements" [4lj . 

III. MODEL FOR THE FMO COMPLEX 

Much of the "quantum biological" interest has fo- 
cused on energy transport in the Fenna-Matthews-Olson 
(FMO) protein complex *23], which is believed to be the 
main contributor to ultra-efficient energy transfer in pho- 
tosynthesis. The FMO protein complex is a trimer in the 
bacterial species prosthecochloris aestuarii. The theoret- 
ical models in the literature I, II E, El, El apply to the 
dynamics of one unit of the trimer. The model assumes 
that a photon impinges on the peripheral antenna of the 
light harvesting complex. Absorption of the photon pro- 
duces an electronic excitation, an exciton, that then tra- 
verses a network of seven chromophores or sites in one 
unit of the trimer. The exciton can either recombinc, 
representing a loss of the excitation, or it can transfer to 
a reaction center, where a light-to-charge conversion oc- 
curs before energy storage. Theoretical models [l2l . [l3[ 
indicate that coherent quantum effects combined with 
decoherence may lead to a quantum stochastic walk [42j 
that transports energy efficiently. Rebentrost et al. pro- 
vide evidence that coherent quantum effects are respon- 
sible for the ultra-high efficiency of photosynthesis [l3[ 
by demonstrating that the transport efficiency is much 
higher with coherent quantum effects than it is without. 

Our physical model for excitation transfer is the nine- 
level model in Refs. [3, m a modification of the seven- 
level model in Refs. [121 fill, Hr3 |. We can restrict dy- 
namics to the single-excitation manifold because the ex- 
citation number is a conserved quantity in the absence 
of light-matter interaction events (within the exciton 
recombination time scale of 1 ns [4J]). The possible 
states for the exciton can be expressed in the site ba- 
sis {\m)y m=1 , where the state |m) indicates that the ex- 
citation is present at site m. The incoherent dynamics 
include a "ground" state \G) corresponding to the loss 
or recombination of the excitation and a sink state \S) 
corresponding to the trapping of the exciton at the reac- 
tion center. The excitation evolves into one of the two 
states \G) or \S) in the limit of infinite time. The den- 
sity operator p for this open quantum system admits the 
following representation in the site basis: 

m,nG{G,l,...,7,S} 

We simplify our analysis by assuming that the dy- 



namics of the density operator are Markovian. Thus, 
a Lindblad master equation, with coherent and incoher- 
ent components, models the dynamics df|[46[. Coherent 
evolution occurs according to the following Hamiltonian 
H 0: 



H = 
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-15.1 
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330 


32.7 


-7.8 
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5.1 


-61.5 


-2.5 


32.7 


280 



(0) 

where the above matrix representation of H is with re- 
spect to the site basis {|w)} m=1 , and the units of energy 
are cm - 1 (the typical units of choice in spectroscopy ex- 
periments). The diagonal terms in H correspond to the 
site energies, and the off-diagonal terms correspond to 
intersite couplings. The order of magnitude of the ener- 
gies in the above Hamiltonian is 100 cm -1 , implying that 
we should observe dynamics on the order of 300 fem- 
toseconds (fs) (Fleming et al. experimentally observed 
behavior on this order in Ref. [171]). 

Three Lindblad superoperators [4^, |46| also contribute 
to the dynamics of the density operator in the nine-level 
model in Refs. [!, H[. The general form of a Lindblad 
superoperator C (p) is 

£ (P) = E Cm (?A mP Al - {AlA m ,p}) , 

m 

where C, m is a rate and A m is a Lindblad operator [45l [it! . 

The first Lindblad superoperator £diss m our model 
corresponds to the dissipative recombination of the exci- 
ton (loss of energy in the system): 

7 

£ diss (p) = £ T m (2 \G) (m| p |m) (G\ - {\m) (m\ , p}) . 

rn — 1 

An excitation at site |m) recombines with rate r m , and 
\G) (m| is the Lindblad operator that effects this dissipa- 
tion. 

The next Lindblad superoperator £ s i n k (p) accounts for 
the trapping of the exciton in the reaction center: 

Ai„k (P) = T sink (2 \S) (3| p |3} (S\ - {|3) (3| , p}) . (8) 

The Lindblad superoperator £ s j n k includes the operator 
\S) (3 1 because evidence suggests that site 3 in the FMO 
complex plays a crucial role in transferring the exciton 
to the reaction center [24j, where it is later exploited for 
energy storage. 

The final Lindblad superoperator £deph accounts for 
the unavoidable dephasing interaction with the environ- 



4 



ment: 



— 1 otherwise. Let Q\m denote the resulting observable 
where 



Aieph (p) = Tm (2 H (m| p |m) (m| - {|m) (m| , p}) , Q w = |^) (^| - (7 - |V>> (^|) = 2 |V> (VI 



7. 



(11) 



( 9 ) 

where 7 m is the rate of dephasing at site m. Discussion of 
this rough treatment of the decoherence appears in Refs. 

The following Lindblad quantum master equation gov- 
erns the evolution of the density operator p: 

P = —i [H, p] + £ d iss (p) + Arink (p) + £dcph (p) , (10) 

where we explicitly see the contribution of the Hamilto- 
nian ([6]) and the noise superoperators ([7][9]) to the dy- 
namics, and we implicitly set h = 1. Evolution according 
to the above Lindblad evolution equation is completely- 
positive and trace-preserving (CPTP) for any time [45j |. 
and we let A/t,t (p) denote the induced CPTP map corre- 
sponding to the evolution of the density operator p from 
an initial time £q to some later time £. 



We might build a dichotomic observable from states in 
the exciton basis. The exciton basis is the energy eigen- 
basis {|^m)}m=i °f the Hamiltonian H in ©, where 



Vm 



E„ 



Then the dichotomic observable constructed from an en- 
ergy eigenstate is Q\<p m )- Note that observables of this 
form commute with the Hamiltonian H in ©■ 

Another possibility is to build the dichotomic observ- 
able from the site basis. This type of observable asks 
the question, Is the excitation at site m ? where m € 
{G, 1, . . . , 7, S}. The dichotomic observable constructed 
from a site state is the site observable Q\ m \- Observables 
of this form do not commute with the Hamiltonian H in 

©■ 



B. Analytical Results for Coherent Dynamics 



IV. THE LEGGETT-GARG INEQUALITY AND 
THE FMO COMPLEX 

In the forthcoming subsections, we consider the appli- 
cation of the Leggett-Garg inequality to the FMO com- 
plex. We first discuss several observables that one might 
measure in a Leggett-Garg protocol. We then obtain an- 
alytical results when the dynamics are purely coherent. 
These analytical results allow us to determine the time 
intervals between measurements in a Leggett-Garg pro- 
tocol that lead to the stongest violation of the inequality. 
We finally use these time intervals in a numerical sim- 
ulation of the FMO dynamics that includes the effects 
of noise. The result is that several observables exhibit 
a strong violation of the inequality for temperatures be- 
low room temperature, and the violation persists in some 
cases up to room temperature. 



A. Observable for the Leggett-Garg Inequality in 
the FMO Complex 

Recall that the Leggett-Garg quantity in ([TJ) involves 
any three dichotomic observables Q(t\), Q{ti), and 
Q (£ 3 ) measured at respective times t\, £2, and £3. 

We have freedom in choosing both the observables that 
we measure and the times at which we measure them. 
Perhaps the simplest dichotomic observable that we can 
construct corresponds to the question |4l| : 

"Is the system in state \ip) or not?" 

The two-element set of measurement operators for this 
question are as follows: {\ip) (ip\ ,1 — We assign 

the value +1 if the system is in the state and the value 



Let us first suppose that evolution of the FMO complex 
is coherent and does not include the noisy Lindblad evo- 
lution terms in (JTHHl). This assumption is unrealistic, but 
it gives a starting point for understanding the Leggett- 
Garg inequality and the FMO complex before proceeding 
with the full-blown evolution in (fTT))) . 

The general form of the two-time correlation functions 
appearing in the Leggett-Garg quantity K in |T]) are as 
follows for a coherent evolution: 



' (t 2 ) e-^-* 1 ) {Q (t0 , p t J e^ 2 " 41 ) 
•{h)e-^-^{Q{t 2 ) lPt2 }e m ^-^ 
'(i 3 )e-^ t3 -* 1 ){g(t 1 ),p tl }e l ^" t2 ) 



C lj2 = iTr 
C 2 , 3 = ilY 

C M = It- 

where H is the Hamiltonian in ([B]), p tl is the ini- 
tial density operator of the FMO complex, p t2 = 
e -»ff(t2-ti)p ti giH(t2-ti) ) an d { v j. j s tne anticommuta- 

tor. As indicated in Section [TTl the correlator Cx,a char- 
acterizes an experiment in which measurements occur at 
times t\ and £3 but no measurement occurs at time £2- 

Suppose first that we prepare the FMO complex in 
some state and take the Leggett-Garg observable to 
be Qu\ as in (llip . Thus, this measurement is a "sur- 
vival probability" measurement (4L|, IH, H^]. Suppose 



further that we measure 



at uniform time intervals 



A£. A straightforward calculation [48|, |49| shows that the 
Leggett-Garg quantity in © is equal to 



4|(V#2At>r - 4 Re («#A t }f <^2AtW 



(12) 



where \ipt) = e~ l \ip). Recall that a violation of the 
Leggett-Garg inequality occurs when the above quantity 
drops below zero. 
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FIG. 1: (Color online) The Leggett-Garg quantity as a function of the uniform time interval At for all seven sites in the FMO 
complex when evolution is coherent and the initial state is the classical mixture 7ri,6- The units of At are picoseconds. A dotted 
line divides each plot into two regions. Points above a dotted line are in a "no violation" region, while points below are in a 
"violation" region. The convention is the same in Figures [2] and [4] To show overall behavior, the At axis has a coarse scale 
up to five picoseconds (ps). However, an experimental test would require control of At to hundreds of femtoseconds. 
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TABLE I: The first column lists the initial state of the FMO complex. The second column lists the site that we measure in a 
Leggett-Garg protocol. The third column gives the analytical form of the corresponding Leggett-Garg quantity as a function 
of the uniform time interval At when dynamics are purely coherent. We use these formulas to compute the results of Figure [1] 
and Table HU 



If the state \ip) is an eigenstate \4> m ) in the exciton 
basis we should not expect to violate the Leggett-Garg 
inequality because the observable Qu m \ commutes with 
the Hamiltonian. We can confirm this intuition by plug- 
ging the eigenstate \<j) m ) into (fT^I) . Doing so gives a value 
of zero for (|12p. thereby saturating the Leggett-Garg in- 
equality, but yielding no violation. 

Thus, we set the Leggett-Garg observable to a site ob- 
servable Q\ m ) where m S {1, ... ,7}. Sites 1 and 6 are 



the chromophores that are closest to the chlorosome an- 
tenna and are thus most likely to be the initial state of 
the FMO complex The initial state p tl of the FMO 
complex can be a pure state |1) (1| or |6) (6| or a uniform 
classical mixture of sites |1) and |6): 

ti, 6 = ^(|1)(1| + |6)<6|). (13) 
One might also consider setting the initial state to the 
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TABLE II: The first column lists the initial state of the FMO 
complex. The second column lists the site observable that the 
Leggett-Garg protocol measures. The third column lists the 
strongest violations for each site observable, and the fourth 
column gives the corresponding time interval At that leads 
to this violation. (For comparison, the strongest possible vi- 
olation of the inequality is —0.5.) We obtained these values 
assuming that evolution is coherent (though, we examined 
times up to At — 5 ps only). The last row in the table indi- 
cates that we do not obtain a violation for any site observable 
when the initial state is the maximally mixed state. 

maximally mixed state 

^ = \(\1) (1I+--- + I7) (7|). 

TableUJlists exact expressions for the Leggett-Garg quan- 
tity in (J3J) for these different cases. 

Figure [T] plots the Leggett-Garg quantity for each site 
observable as a function of the uniform time interval At, 
when the initial state is 711,6 • The result is that each 
of the seven site measurements gives a violation of the 
Leggett-Garg inequality for some intervals At when the 
dynamics are purely coherent. These results may not be 
particularly surprising [48l 143 ] given that the system is 
quantum, the measurements are sharp, and the dynamics 
are purely coherent. 

Table [TT] lists the time intervals that lead to the 
strongest violation for At £ [0, 5] ps. We only consider 
times up to five picoseconds (ps) because it is likely that 



the exciton will trap by this time for the case of inco- 
herent dynamics in the next section. Table [TT] also lists 
the time intervals that lead to the strongest violation 
when the initial state is |1) (1|, |6) (6|, or tt. We now 
use these time intervals for the more realistic case, where 
the excitation in the FMO complex experiences the noisy 
contributions in (J7J15]) from dissipation and dephasing. 

C. Numerical Results for Incoherent Dynamics 

The Lindblad evolution in (flQ|) is a more realistic 
model for evolution of the excitation in the FMO com- 
plex. It incorporates the effects of excitonic recombina- 
tion, trapping to the reaction center, and enviromental 
dephasing noise, albeit in a Markovian context. The two- 
time correlation functions appearing in the Leggett-Garg 
quantity K in (JTJ) must now be evaluated using the fol- 
lowing more general forms: 

Ci, 2 = ~Tr[Q (t 2 )Af t2M ({Q(*i), />})], 

C 2 ,3 = \tt [Q {{Q (t 2 ) ,M 2 , tl (P)})] , 

C 1>3 = iTr[Q (t 3 )Af t3 , tl ({Q(ti), /»})], 

where Af is the superoperator that propagates the density 
operator forward in time according to the evolution in 

In order to perform a numerical simulation for the in- 
coherent dynamics (|10p . we need to fix the trapping rate 
r s j n k, the recombination rates T m , and the dephasing 
rates 7™. We take these rates from Refs. [3, 5|. Un- 
less otherwise stated, the trapping rate = 62.8/1.88 
cm -1 , corresponding to about 2ir ■ c ■ T s i n ^ w 6 ps -1 , 
where c is the speed of light in units of cm-ps -1 . The 
recombination rates T m are uniform for all sites in the 
FMO complex so that T m = T lecomh = 1/ (2 • 188) cm" 1 , 
corresponding to about 2ir ■ c ■ r rec omb ~ 5 x 10~ 4 ps -1 . 
We assume that the dephasing rate is uniform so that 
7m = 7 for all sites m. 

As in the previous subsection, we choose the Leggett- 
Garg observable to be a site observable Q\ rn ) where 
m £ {1, . . . , 7}. The initial state pt 1 can either be the 
pure state |1) (1|, |6) (6|, or a uniform classical mixture 
of sites |1) and |6) as in (jT3J) , for the same reasons men- 
tioned in the previous subsection (we do not consider the 
maximally mixed state because there is no phenomeno- 
logical evidence for this case, and furthermore, it does 
not give a violation even in the coherent case) . The time 
interval At between measurements of the Leggett-Garg 
observable is taken from Table HT1 FigureslHO andHJdis- 
play the Leggett-Garg quantity as a function of the de- 
phasing parameter 7 for each site observable Q\ rn ), where 
m£ {1,...,7}. 

The figures demonstrate that several site observables 
exhibit a violation even as 7 increases. The temperature 
analysis of the environment in Ref. [l3j indicates that a 
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Dephasing rate y in units of inverse picoseconds 

FIG. 2: (Color online) The above figure displays the Leggett- 
Garg quantity as a function of the dephasing parameter 7 
for each site observable Q\ m ), where m G {1, . . . , 7}. The 
site observables Q\ 3 \ and Qu) do not give a violation of the 
Leggett-Garg inequality for any amount of dephasing noise, 
even though we observed a violation in Figure [1] with purely 
coherent dynamics. Site observables Qu), Q|2)> Q|5)> Q|6)i 
and Q| 7 ) give a strong violation of the inequality for low 
dephasing noise. Each of their corresponding Leggett-Garg 
quantities has a smooth, monotonic transition to the "no vio- 
lation" region for stronger dephasing noise, with the Leggett- 
Garg quantity for site observables Q|i), Q\2), Q\s), Q\e) with- 
standing the strongest amount of dephasing noise before they 
make a transition to the "no violation" region. Using the tem- 
perature analysis of the environment in Ref. [l3T |. a dephasing 
rate of 2.1 ps _1 corresponds to a temperature of around 77°K 
and a dephasing rate of 9.1 ps _1 corresponds to a tempera- 
ture of 298°K (room temperature). Site observables Q\i), 
Q\2), Q\5), an d Q|6> give a violation for room temperature, 
with the Leggett-Garg quantity respectively equal to —0.0039, 
-0.0015, -0.0059, and -0.0079. It might be difficult for an 
experimentalist to observe these violations at room tempera- 
ture given that they are weak. 




Dephasing rate y in units of inverse picoseconds 



FIG. 3: (Color online) This figure is similar to Figure [2] with 
the exception that the initial state of the FMO complex is a 
pure state at site 1. Site observables Qm\ and Q| 2 > give a vio- 
lation for room temperature, with the Leggett-Garg quantity 
respectively equal to —0.0077 and —0.003. 




Dephasing rate y in units of inverse picoseconds 



dephasing rate of 2.1 ps _1 corresponds to a temperature 
of around 77°K and a dephasing rate of 9.1 ps _1 corre- 
sponds to a temperature of 298°K (room temperature). 
Under the assumption that this is approximately correct, 
Figures [H [3J and Q] predict a violation at both tempera- 
tures. However, it might be difficult for an experimental- 
ist to observe these violations at room temperature given 
that they are weak. 

We have verified the robustness of these predictions 
against variations in the dynamical parameters in (|10[) . 
Ref. [24| mentions that the site energies of their calcu- 
lated Hamiltonian are accurate within ±2 cm -1 . We 
therefore conducted several simulations that added in- 
dependent, zero- mean Gaussian noise with variance 2 to 
each entry in the Hamiltonian matrix in ^ to deter- 



FIG. 4: (Color online) This figure is similar to Figure [2] with 
the exception that the initial state of the FMO complex is a 
pure state at site 6. Site observable gives a violation for 
room temperature, with the Leggett-Garg quantity equal to 
-0.0155. 



mine if the violations would still hold under this slight 
perturbation. The result is that all room temperature vi- 
olations in Figures [21 [3j and [4] still hold, and in fact, the 
values of the Leggett-Garg quantities are the same up to 
the fourth decimal place. The trapping rates in the liter- 
ature vary substantially, including 0.25 ps _1 [24[, 1 ps" 1 
[13], or 4 ps _1 16], so we have checked our results for all 
of these choices. The plots for all these cases are similar 
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to Figures [3l and [4] and have approximately the same 
values for violations at room temperature. This finding 
is reasonable given that most of our measurement times 
are less than the average trapping times corresponding 
to these other trapping rates. 

Our numerical simulations demonstrate that several 
choices of measurements lead to a violation of the 
Leggett-Garg inequality even when noise processes act on 
an excitation in the FMO complex. Experimental confir- 
mation could irrevocably exclude a class of macrorealistic 
theories from describing the dynamics of the excitation. 

V. CONCLUSIONS 

We have framed tests that could be used to experimen- 
tally exclude a class of macrorealistic theories, including 
a classical incoherent hoppping model, from describing 
the room-temperature dynamical behavior of an excita- 
tion in the FMO complex. To do so, we have introduced 
several examples of observables that one might apply in 
a test of macrorealism, and our numerical simulations 
predict that these observables lead to a violation of the 
Leggett-Garg inequality. 

We discuss several ideas for furthering this research. 
The Leggett-Garg measurements that we have consid- 
ered here are projective measurements. In practice, an 
experimentalist is unlikely to realize such an idealized 
measurement. An experimentalist is more likely able 
only to effect a "noisy" or "fuzzy" measurement of the 
sites in the FMO complex. Along these lines, Kofler et 
al. have shown that "fuzzy" measurements may lead 
to our observation of a classical world even in the pres- 
ence of quantum dynamics [4^, [4^ |. Future work should 



formulate Leggett-Garg inequalities involving quantities 
that are conveniently accessible to FMO experimentalists 
(measurements in the site basis are currently infeasible 
|50() and should consider the limits of experimental mea- 
surement capabilities. Furthermore, future work should 
consider more realistic models of noise in the FMO com- 
plex, potentially inclu ding correlated noise (25l . l5l| and 
non-Markovian effects [26| . 

One might consider an application of Bell's inequal- 
ity to study the ability of the FMO complex to preserve 
entanglement. 

We have applied the Leggett-Garg inequality to the 
FMO complex, and it would be valuable to apply the 
Leggett-Garg test to the dynamics of magnetoreception 
in the avian compass [1, d, 0, [l(| • It would also be inter- 
esting to explore the Leggett-Garg inequality in artificial 
quantum networks with particular noisy interactions. A 
study of this sort might lead to an increased understand- 
ing of the dividing line between macrorealism and non- 
classicality in more complicated systems. One might be 
able to apply the ideas in Refs. [H, El El, El] here. 
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